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ABSTRACT 

We present some results of a new calculation of the 0(a) 
electroweak radiative corrections to W boson production at 
hadron colliders with special emphasis on the transverse mass 
distribution. 

I. INTRODUCTION 

Despite the remarkable success of the Minimal Standard 
Model (MSM) in describing elementary particle interactions at 
presently accessible energies, there is little direct experimental 
information on the mechanism which generates mass for the W 
and Z bosons. In the MSM, spontaneous electroweak symme- 
try breaking is responsible for mass generation. The existence 
of the Higgs boson in the MSM is a direct consequence of this 
mechanism. 

Complementary to the direct Higgs boson search at colliders 
(Mjj > 58.4 GeV JjJ] from LEP1, where Mjj is the Higgs bo- 
son mass), indirect information on Mjj can be extracted by con- 
fronting theoretical predictions for radiative corrections to elec- 
troweak observables with high precision measurements. As- 
suming that the MSM is valid, a global fit to the currently avail- 
able data from LEP and SLC with a s , the top quark mass (m t ) 
and Mjj as free parameters yields Mjj = 146+gg 2 GeV [§]. 
Similar results have been obtained in Ref. [^J. The indirect con- 
straints on Mjj are expected to improve considerably in the fu- 
ture with more precise measurements of the top quark mass and 
the W boson mass (My/)- Presently, their world averages are 
m t = 175 ± 6 GeV g] and M w = 80.356 ± 0.125 GeV [§. 
The precise measurement of Mw is therefore one of the priori- 
ties of future collider experiments. LEP2 and Runll at Fermilab 
(JCdt = 2 fb _1 ) are aiming for an uncertainty on Mw of about 
40 MeV [||] and 35 MeV (per experiment) [0], respectively. 
Further upgrades of the Tevatron accelerator complex (TeV33) 
could yield an overall integrated luminosity of O(30 fb _1 ), and 
a precision of Mw of about 15 MeV m[. Finally, it may be pos- 
sible to measure Mw at the LHC with an accuracy better than 
15 MeV, if an integrated luminosity of 10 fb -1 is accumulated 
with the accelerator operating at a reduced luminosity of about 
10 33 cm" 2 s- 1 (see Ref. [§]). 

Obviously, to measure Mw with high precision, it is crucial to 
fully control higher order electroweak (EW) and QCD radiative 
corrections. 

In this contribution we present some results of a new cal- 
culation of the EW 0(a) corrections to W boson production 
in hadronic collisions. In particular, we study the effect of 
these corrections on the W transverse mass (My(fo')) distri- 



bution from which Mw is extracted at the Tevatron. In a pre- 
vious calculation Jic|], only the final state photonic corrections 
had been included, using an approximation in which the sum 
of the virtual and soft part is indirectly estimated from the 
inclusive 0(a 2 ) W — » tv(^[) width and the final state hard 
bremsstrahlung contribution. Our calculation includes both ini- 
tial and (complete) final state corrections, as well as their inter- 
ference. As a result of our calculation, the current systematic 
uncertainty of AM W = 20 MeV [[yj ^ originating from EW 
radiative corrections will be reduced. We shall only discuss the 
W — > ev e decay channel here. More details and a discussion of 
the W — > [iv ^ decay channel will be presented elsewhere [pj||. 

II. THE 0(a) CONTRIBUTION TO 
RESONANT W PRODUCTION 

When calculating the EW radiative corrections to resonant W 
boson production the problem arises how to treat an unstable 
charged gauge boson consistently in the framework of perturba- 
tion theory. This problem has been studied in Ref. [ 14 1 with par- 
ticular emphasis on finding a gauge invariant decomposition of 
the EW 0(a) corrections into a QED-like and a modified weak 
part. Unlike the Z boson case, the Feynman diagrams which 
involve a virtual photon do not represent a gauge invariant sub- 
set here. In Ref. [|4D, it was demonstrated that gauge invariant 
contributions can be extracted from the infrared (IR) singular 
virtual photonic corrections. These contributions can be com- 
bined with the real photon corrections in the soft photon region 
to form gauge invariant QED-like contributions corresponding 
to initial state, final state and interference corrections. The soft 
photon region is defined by requiring that the photon energy in 
the parton center of mass frame, E 1 , is smaller than a cutoff 
AE = 5 s V§/2, where is the parton center of mass energy. 
In this phase space region, the soft photon approximation can 
be used to calculate the cross section as long as S s is sufficiently 
small. Throughout the calculation the soft singularities have 
been regularized by giving the photon a fictitious mass. In the 
sum of the virtual and soft photon terms the unphysical photon 
mass dependence cancels, and the QED-like contributions are 
IR finite. 

The IR finite remainder of the virtual photonic corrections 
and the weak one-loop corrections can be combined to sepa- 
rately gauge invariant modified weak contributions to the W 
boson production and decay process. Both the QED-like and 
the modified weak contributions can be expressed in terms of 
form factors, Fq ED and F® eak , which multiply the Born cross 



section [n_4J. The superscript a in the form factors denotes the 
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initial state, final state or interference contributions. 

The complete 0(a 3 ) parton level cross section of resonant 
W production via the Drell-Yan mechanism q{q v — > ff'i'j) 
can then be written as follows [|l4]] : 



[d&^ F* ED (s,t) + da^ 3 ] , (1) 
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where the Born cross section, da^°\ is of Breit-Wigner form 
and s and t are the usual Mandelstam variables in the parton 
center of mass frame. The modified weak contributions have to 
be evaluated at s = Jl4| ] . The IR finite contribution da^^ 
describes real photon radiation with E 1 > AE. 

Additional singularities occur when the photon is emitted 
collinear with one of the charged fermions. These collinear 
singularities have been regularized by retaining finite fermion 
masses. Thus, both da^^ z and Fq ED (a — initial, final) 
contain large mass singular logarithms which have to be treated 
with special care. In the case of final state photon radiation, the 
mass singular logarithms cancel when inclusive observables are 
considered (KLN theorem). For exclusive quantities, however, 
these logarithms can result in large corrections, and it may be 
necessary to perform a resummation of the soft and/or collinear 
photon emission terms. For initial state photonic corrections, 
the mass singular logarithms survive. These terms are universal 
to all orders in perturbation theory, and can therefore be can- 
celled by universal collinear counterterms generated by 'renor- 
malizing' the parton distribution functions (PDF), in complete 
analogy to gluon emission in QCD. 

To increase the numerical stability, it is advantageous to ex- 
tract the collinear part from d&%^ 3 for both the initial and final 
state, and perform the cancellation of the mass singular loga- 
rithms analytically. The collinear region is defined by requiring 
that the angle between the fermion and the emitted photon is 
smaller than a cutoff parameter Sg. The reduced 2^3 con- 
tribution, which comprises the real photon contribution away 
from the soft and collinear region, is evaluated numerically us- 
ing standard Monte Carlo techniques. Our method is very simi- 
lar to the phase space slicing method described in Ref . fll5| ] . 

It should be noted that the analytic cancellation of the final 
state mass singular logarithms is only possible when realistic 
experimental electron identification requirements are taken into 
account. This will be discussed in more detail in Sec. II. B. 

In the remainder of this section, we extract the collinear be- 
havior of d&2^ 3 for both the initial and final state, and perform 
the 'renormalization' of the PDF. 



A. Initial state photon radiation 

In the collinear region, for sufficiently small values of Sg, the 
leading pole approximation can be used, and the initial state real 
photon contribution can be written as a convolution of the Born 
cross section with a collinear factor (see also Ref. [[160): 
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i.e. an incoming parton qi with momentum pi, mass rrii and 
charge quantum number Qi splits into a parton q^ with momen- 
tum zpi and a photon of momentum (1 — z)pi. The value for 
the upper limit of the z integral avoids overlapping with the soft 
photon region. 

The counterterms generated by the 'renormalization' of the 
PDF can simply be derived from Ref. [|l7| by performing the 
replacement 
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One then finds: 
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fv+s = 9 + ^- + 3 log S s - 2 log 2 S s (4) 
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The qi(x) and qt(x, Q 2 ) are the unrenormalized and renormal- 
ized parton distribution functions, respectively. The parame- 
ter Xpc distinguishes between the MS (Xfc — 0) and the DIS 
scheme (Xfc = 1)- The scheme dependent contributions f c and 
f v +s can be derived from Ref. jl8|]. Q is the factorization scale. 

The cross section foxpp — > W{^) — * is then obtained 

in two steps: first, the parton level cross section of Eq. ([[]) is con- 
voluted with the unrenormalized PDF qi(x), and second, qi(x) 
is replaced by the renormalized PDF (x, Q 2 ) by using Eq. (|j). 
The initial state QED-like contribution F'Q E l E al and the collinear 
part da' l ™j t l lal , including the effect of mass factorization, can be 
grouped into a single 2-^2 contribution: 
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As expected, the mass singular logarithms cancel completely. 

In our calculation, we have not taken into account the QED 
radiative corrections to the Gripov-Lipatov-Altarelli-Parisi evo- 
lution of the PDF. This introduces an uncertainty that needs to 
be quantified. We will address this question in Ref. [Oh. A 
study of the effect of QED on the evolution indicates that the 
change in the scale dependence of the PDF is small [H. To 
treat the QED radiative corrections in a consistent way, they 
should be incorporated in the global fitting of the PDF. 

B. Final state photon radiation 

Similar to initial state radiation, the final state real photon cor- 
rections in the collinear region can be described as a convolu- 
tion of the Born cross section with a collinear factor. Using the 
leading pole approximation one finds (see also Ref. [|l6|]): 
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When realistic experimental conditions are taken into account 
(see Sec. Ill), the electron and photon four-momentum vectors 
are recombined to an effective electron four-momentum vector 
if their separation Ai? e7 in the azimuthal angle - pseudorapidity 
plane is smaller than a critical value R c . If the cutoff parameter 
Sg is chosen to be smaller than R c the integration over the mo- 
mentum fraction z in Eq. (0) can then be performed analytically. 
In this case, the mass singular logarithms in the sum of do co n 
and the QED-like contribution Fq££ explicitly cancel, and one 
obtains: 

dolZt = ^2 J dx x dx 2 [qi{xi,Q 2 ) q l ,(x 2 ,Q 2 ) + (1 <-> 2)] 
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The approximation [|10p used so far in modeling the EW radia- 
tive corrections to W boson production at the Tevatron differs 
from our calculation only in the 2^2 contribution. At the par- 



ton level, the difference between Eq. 
is given by 



and the approximation 
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As we shall see, this difference has a non-negligible effect on 
the transverse mass distribution, and thus on the W boson mass 
extracted from experiment. 

m. NUMERICAL IMPACT ON THE 
TRANSVERSE MASS DISTRIBUTION 

Since detectors at the Tevatron cannot directly detect the neu- 
trino produced in the leptonic W boson decay, W — > ev e , and 
cannot measure the longitudinal component of the recoil mo- 
mentum, there is insufficient information to reconstruct the in- 
variant mass of the W boson. Instead, the transverse mass dis- 
tribution of the final state lepton pair is used to extract Mw- 
In the following, we therefore focus on the Mr(et'e) distribu- 
tion. The following acceptance cuts and electron identification 
criteria are used to simulate detector response: 

• the uncertainty in the energy measurement is simulated by 
Gaussian smearing of the lepton momenta corresponding 
to the D0 electromagnetic energy and missing transverse 
energy resolution, 

• the photon and electron are treated as separate particles 
only if Ai? e7 > 0.3. For smaller values of Ai? e7 , the 
four momentum vectors of the two particles are combined 
to an effective electron momentum four-vector (recombi- 
nation cut). In the region 0.2 < AR ei < 0.4 the event is 
rejected if E 1 > 0.15E e (isolation cut), 

• we impose a cut on the electron transverse energy of 25 
GeV, a missing transverse energy cut of 25 GeV, and re- 
quire the electron pseudorapidity to be \rj e \ < 1.2. 

We use the MRSA set of partem distributions [^0|] with the fac- 
torization and renormalization scales set equal to Mw- For 
the numerical evaluation of da^fi 0,1 ( see Eq. @) we use the 
MS scheme. The leading order (LO) and next-to-leading order 
(NLO) Mx{ev) differential cross sections are shown in Fig. |l| 
As can be seen, the overall effect of the 0(a) corrections is to 
reduce the cross section. 

The dependence on the cutoff parameters 6 S and 6$ must can- 
cel in the sum of the 2 — > 2 and reduced 2 — > 3 contributions, 
provided that the cutoff parameters are chosen sufficiently small 
such that the soft photon and leading-pole approximation are 
valid. This is demonstrated in Fig. || for the final state contribu- 
tions. We have also verified the cancellation for the initial state 
and interference contributions. 

In Fig. ||, the effect of the various individual contributions to 
the EW 0(a) corrections on the Mr(ei/ e ) distribution is shown. 
The following conclusions can be drawn: 
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Figure 1: The Mx{ev e ) distribution of W boson production at 
the Tevatron (y/S = 1.8 TeV). 
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• The initial state QED-like contribution uniformly increases 
the cross section by about 1%, but is largely cancelled by 
the modified weak initial state contribution. 

• Both the complete O(a) initial state contribution and the 
interference contribution are small and change the shape 
of the Mt(&v) distribution very little. These contributions 
therefore, are expected to have a negligible effect on the 
value of Mw extracted from the data. 

• The final state QED-like contribution significantly changes 
the shape of the transverse mass distribution and will, 
therefore, have a non-negligible effect on the value of Mw 
extracted from data. As for the initial state, the modified 
weak final state contribution reduces the cross section by 
about 1%, but has only a small effect on the shape of the 
transverse mass distribution. 

The change in the shape of the Mr(e^ e ) distribution due to the 
QED-like final state corrections can be easily understood. Pho- 
ton radiation reduces the energy of the final state electron and 
thus the transverse mass when the electron and photon momenta 
are not recombined. 

In Fig. we display the ratio of the Mr(ez/ e ) distribution 
obtained with our complete NLO calculation to the one ob- 
tained by using the approximation of Ref. [|To|]. The depen- 
dence of the ratio on the transverse mass is described by Eq. (Q). 
For Mx{ev) < Mw, most events originate from the region 
s « M w , due to the Breit-Wigner resonance. Consequently, 
there is very little dependence on Mt in that region. For 
Mt(zv) > M\y, the steeply falling cross section in the tail of 
the Breit-Wigner resonance favors events with s w Mf.. There- 
fore, in this region, the terms proportional to \og(s/M w ) in 
Eq. (H) cause a change in the shape of the transverse mass dis- 



Figure 2: The 2-^2 and reduced 2 — > 3 final state contribution 
a) as a function of the cutoff S s with 8g = 0.001, and b) as a 
function of Sg with S s = 0.01. 



tribution. The difference in the line shape occurs in a region of 
the Mr(ei/ e ) distribution which is sensitive to Mw, and we ex- 
pect that the W boson mass extracted when using the complete 
calculation will differ by several MeV from the value obtained 
using the approximate calculation. Since this difference is much 
smaller than the present uncertainty for Mw, the approximation 
of Ref. [ |lO| ] provides an adequate description of W boson pro- 
duction at the Tevatron for the currently available data. How- 
ever, for future precision experiments, a difference of several 
MeV in the extracted value of Mw can no longer be ignored, 
and the complete calculation should be used. 

IV. CONCLUSIONS 

We have presented the results of a calculation of the EWO(a) 
corrections to W boson production at hadron colliders. Both 
initial and (complete) final state corrections, as well as the in- 
terference between the initial and final state corrections are in- 
cluded in our calculation. The initial state corrections and the 
interference contribution are found to be small and uniform in 
Air(ei'e), and are expected to have a small effect on the W bo- 
son mass extracted from experimental data. As expected, the 
final state corrections dominate the EW radiative corrections. 
They significantly change the shape of the transverse mass dis- 
tribution, and thus the value of Mw extracted from the data. 

We also compared the result of our complete calculation with 
that of Ref. [ 1 1 which uses an approximation to estimate the 
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Figure 3: The ratio of the NLO to LO Mr(ei/ e ) distribution for 
various individual contributions: the QED-like initial or final 
state contributions (solid), the complete 0(a) initial and final 
state contributions (short dashed) and the initial-final state in- 
terference contribution (long dashed). 
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Figure 4: The ratio of the complete to the approximate NLO 
MT{ev e ) distribution. 



sum of soft and virtual final state corrections. For Mt > My/, 
the complete and approximate NLO A/r(e^ e ) distributions dif- 
fer in shape, and we expect that they will yield values for the W 
boson mass which differ by several MeV. 

Our calculation substantially improves the treatment of EW 
radiative corrections to W boson production in hadronic colli- 



sions, and will allow to significantly reduce the systematic un- 
certainties associated with these corrections when the W boson 
mass is extracted from Tevatron data. 
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